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1.  Introduction 


T 


From  potential  theory  it  is  well  known  that  a series  of  spherical 
harmonics  of  the  gravity  field  is  convergent  outside  the  "minimum  sphere" 
enveloping  all  masses.  In  many  studies  it  has  been  shown  that  for  certain 
earth  models  the  series  may  be  extended  to  a sphere  entirely  within  the  surface 
of  matter  or  the  radius  of  convergence  might  be  located  somewhere  between 
these  extremes.  For  references,  see  Moritz  (1961),  Molodensky  et  al.  (1962, 
pp.  118-120),  Pick  (1965),  Morrison  (1969),  Hotine  (1969,  pp.  172-173)  and 
Levallois  (1972). 

_ . < 

Moritz  (ibid.)  showed  that  if  the  earth  were  a homogeneous  oblate 
spheroid  ( a level  ellipsoid),  then  the  series  for  V would  converge  at  the  surface. 
Levallois  studied  the  condition  for  extension  of  the  series  down  to  the  surface 
of  a homogeneous,  approximately  spherical  body.  For  such  a body  he  found  the 
following  condition  for  convergence  at  the  surface: 


where 


kB<  (0.132)" 

y (Cj+Sn.^) 


2n+  1 


CflB  > Sna  = fully  normalized  spherical  harmonics 


For  the  eairth  we  have,  according  to  Kaula's  rule*. 


(1.1) 


10"  Vn^ 


Thus  the  condition  above  is  not  satisfied  in  this  case.  Some  other  models  and  the 
corresponding  radii  of  convergence  of  the  series  are  given  by  Morrison  (ibid). 

However,  in  all  these  studies  the  models  are  either  homogeneous  or 
bodies  of  revolution.  Already  a small  disturbing  body  with  mass  centre  located 
outside  a homogeneous  sphere  makes  the  harmonic  series  divergent  at  the 
surfoce  of  the  sphere  (see  the  previous  references  of  Moritz,  Molodensky  et.al. , 
Pick  and  the  example  3.1  below).  Hence,  because  of  the  irregular  mass 
distributions  of  the  actual  earth  the  series  of  spherical  harmonics  must  be 
considered  divergent  at  the  surface  of  the  earth  (Moritz,  1961). 
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These  theoretical  aspects  do  not  imply  that  the  anal3d:ic  continuation 
of  gravimetric  quantities  down  to  the  surfoce  of  the  earth  is  meaningless.  Moritz 
(1969)  has  paid  attention  to  this  fact  and  states  that  asymptotic  series,  which 
are  mathematically  speaking  divergent,  are  frequently  used  in  mathematical 
physics.  "Such  series  can  be  used  if  the  first  terms  decrease  rapidly  enough 
for  their  sum  to  provide  a good  estimation  to  the  function  to  be  calculated;  it 
will  not  matter  practically  if  the  ne^ected  higher  terms  will  start  increasing 
again.  The  practical  use  of  divergent  series  needs,  however,  to  be  justified. 

An  arbitrarily  accurate  approximation  can  be  obtained  only  with  convergent  series; 
with  divergent  asymptotic  series,  this  error  cannot  be  reduced  below  a certain 
limit.  It  must  be  investigated  whether  this  limit  is  small  enough  so  as  to  be  in 
keeping  with  the  desired  accuracy.  " 

The  main  object  of  this  report  is  to  estimate  the  downward  continuation 
error  of  the  representation  of  the  gravity  field  by  a series  of  spherical  harmonics 
at  the  surface  of  the  earth.  We  start  with  a definition  of  the  problem. 


2.  Definition  of  the  Problem 

In  the  volume  external  to  a sphere  enclosing  all  mass  of  the  earth ^ 
(the  minimum  sphere,  or  Brillouin  sphere)  the  gravity  potential  'of  the  earth  is 
harmonic  and  can  be  expanded  into  a series  of  spherical  hairmonics: 


(2.1) 

where 


Vn  = Y (C„,  cosmX  + S„,  sinmX)  P„,(cose) 

R — < 

■ = 0 

G = Newton's  constant  of  gravitation 
M = mass  of  the  earth 
R = radius  of  the  "minimum  sphere" 

(r,  6,X)  = spherical  coordinates 

Cn»  > S„.  = fully  normalized  spherical  harmonic  coefficients 
Pni  (cos  6)  = associated  Legendre  function 

In  many  cases  the  potential  V is  substituted  by  the  disturbing  potential: 


i 


i 


T = V - U 


^ The  atmosphere  of  the  earth  is  not  considered 


i 

j 
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Uo=  Vo  t ^20  U40  )•  In  this  case  the  series  expansion  in  formula  (2. 1) 

starts  at  n = 2 and  the  coefficients  for  Tjo  and  T40  are  now  considered  as 


corrections  to  the  normal  potential. 


Inserting  the  series  expansion  of  T into  the  spherical  approximation 
of  ihe  boundary  condition  of  physical  geodesy  (see  Heiskanen  and  Moritz,  1967, 
p.  88): 


(2.2)  ^ =-|I.  - IX 

or  r 


the  gravity  anomalies  {^)  are  obtained  in  a series  of  the  potential  harmonics: 
(2.3)  Ag  = y 

n=a 

where 

n 

(n-l)y  (Cm  cos  m A + Sn,  sinmA)  Pj,  (cos  6) 

' r\  — » 

■ = 0 


From  (2.2)  it  is  easily  shown  that  r^g  is  harmonic  in  the  same  domain  as  T, 
i.  e.  outside  the  surface  of  the  earth  (if  the  influence  of  the  earth's  atmosphere 
is  neglected). 

The  coefficients  C„,  and  S„,  of  the  series  (2.1)  and  (2.2)  have  been 
determined  to  various  degrees  (N)  from  terrestrial  gravity  observations,  satellite 
observations  and  from  combinations  thereof.  From  these  coefficients  the  external 
gravity  field  can  be  determined  by  a truncated  series: 


(2.4) 


ns:  0 


These  truncated  series  are  theoretically  correct  outside  the  "minimum 
sphere"  (R  = 6384.403  km  for  a = 6378. 140  km,  see  Appendix).  Now  the 
question  arises  whether  the  series  V and  Ag  can  be  analytically  continued  down 
to  the  surface  of  the  earth.  More  precisely:  What  is  the  error  of  such  a repre- 
sentation? Is  it  possible  to  find  simple  correction  terms  ? 
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3.  A Simple  Model 


Before  going  deeper  into  the  problem  of  estimating  the  errors  of  using 
spherical  harmonic  series  at  the  surface  of  the  earth,  we  would  like  to  present 
a simple,  illustrative  model. 


A disturbing  point  mass  m is  located  outside  a homogeneous  sphere 
M with  radius  ro  . The  distance  to  m from  the  center  of  M is  R.  See  Figure 
3.1.  At  a point  P of  distance  I from  m the  disturbing  potential  is: 


where 


Gm 

Newton's  constant  of  gravitation 

the  angle  between  the  radius  vectors  r and  R 

distance  from  the  center  of  M to  P 


Figure  3. 1.  A homogeneous  sphere  with  an  exterior  disturbing 
point  mass  m. 


1 i 

f ! 
! I 


At  a point  outside  a sphere  with  the  same  center  as  M and  radius  R the  closed 
formula  for  T can  be  expanded  in  the  following  series  of  Legendre's  polynomials: 


(3.1) 


T.=  Y (COS0)  , r> 


R 


n=  0 


Applying  this  formula  at  the  surface  of  M we  obtain: 
(3.1’)  T = M.  Y J*'p„(cosi(>) 


1=  o 


This  series  is  divergent  for  N = ® , because  R > ro  . The  correct  value  for  the 
disturbing  potential  at  M is  given  by: 


(3.2) 


* (ro +R^-2ro  R cos  (|) ) 


n = 0 


Thus  we  obtain  the  following  error  of  T : 


€r  (N)  = T-  Ti  = 6T(N)  + er  (N) 


where  6T(N)  is  the  error  of  analytic  continuation  and  eT(N)^is  the  truncation 


error.  The  truncation  error  is  caused  by  the  truncation  of  T at  degree  N: 


(•^J^'p„(cos(p) 


n=N+l 


The  downward  continuation  error  is  defined  by: 

N 

6T(N)  = Y 6T„ 


n=  0 


11 
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where  6T„  is  the  error  caused  by  the  improper  downward  continuation  of  (T.) 
to  the  surface  of  the  main  sphere.  We  have; 


■■i 


where 


and 


— (Te)n“  (Tj)n—  Ca  (Tj)b 

(Ti)„=  ^ (-^J^'PaCCOStt) 


4 


The  coefficient  Cn  is  the  relative  error  of  (Tj)n . Inserting  R = r©  + h we  obtain: 

-l^(2n+l)  f 
' ro ' To 


In  Figure  3.2  the  true  value  of  the  geoidal  heights  (Ti/y)  and  some  estimates  for 
downward  continuation  are  shown,  hi  Figure  3. 3 the  relative  error  Cn  is  illustrated. 
Finally,  in  Figure  3.4  the  relative  error 

N 

6T(N)/y  (TOa 

n=  0 


is  given.  The  figures  show  that  the  relative  errors  are  increasing  with  the 
elevation  (h)  of  the  disturbing  mass  and  with  the  degree  (n  respective  N). 


Figure  3. 3.  Percentage  error  in  the  disturbing  potential  harmonic  of  degree  n caused 

by  the  improper  downward  continuation  to  a sphere  with  a point  mass  located 


Figure  3.4.  Percratage  error  in  the  disturbing  potential  harmonic  expansion  to  degree  N 
caused  by  the  improper  downward  continuation  to  a sphere  with  a point  mass 
located  at  a height  h above  that  sphere  (ro  = 6370  km,  |l)  = 3’). 


i 
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For  r < R we  have  from  (3.2): 


(3.4) 


r < R 


a = 0 


Both  of  ttteee  aeries  are  divergent  for  r = R.  to  this  case  we  obtain  the 
derivative  directly  from  the  closed  formula  (3.2): 


(3.5) 


/9T\  M r = ti 

\SrA=a  2/2  R®(l-co8 


r = R 


The  gravity  anomaly  is  given  by  formula  (2. 2): 


dT  2T 

J7  ~ T 


Hence  the  formulae  for  Ag  corresponding  to  (3.3)  throu^  (3*5)  are: 


(3.6) 


n ^ 8 


R 


ns  0 


(3.7) 


and 


Agi  = -|s  Y (n  + 2)(~-)  ^P,(cos0)  , r<R 


n=  O 


(3.8) 


3/2 


r=  R 


,A_.  U 3''  is  A A, 

(^)r=a  = “ ^ 1 t ^ - 

« 4(l-cos^)e  0 


Let  us  assume  that  we  use  a truncated  form  of  formula  (3. 6)  to  represent  the 
gravity  anomaly  at  the  surtoce  of  the  main  sphere  (r  = ro): 


nsO 


! 

I 
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where  N is  the  degree  of  truncation.  The  true  anomaly  at  this  level  is  given  by 
(3.7).  Thus  the  total  error  of  ^ becomes: 


where  6 ^ (N)  is  the  error  of  the  analytic  continuation  and  e^  (N)  is  the  truncation 
error.  We  obtain: 


where 


The  coefficients  Cb  and  | dn  | are  the  relative  errors  of  the  n th  harmonic  of 
the  anomalous  potential  and  the  gravity  anomaly,  respectively.  Cg  {Approaches 
zero  for  R-^ro . This  is  not  the  case  for  da . The  series  is  not  conver- 
gent for  ro  = R.  The  reader  should  notice  that  the  usual  relation  between  the 
harmonics  of  the  potential  and  the  gravity  anomaly  is  not  valid  in  this  case 
(see  Heiskanenand  Moritz,  1967,  p.  97): 


The  influence  of  the  disturbing  masses  on  the  anomalies  is  more  pronounced 
than  for  the  potentials.  One  reason  for  this  is  that  while  the  disturbing 
potential  is  dependent  only  on  the  distance  to  the  disturbing  masses,  the 
contribution  to  the  gravity  anomalies  from  masses  located  above  the  point 
of  computation  usually  have  opposite  sign  to  that  obtained  in  ihe  downward 
ctmtinuaticm  procedure. 

Finally,  we  derive  the  errors  of  the  vertical  gradient  of  the  gravity 
anomalies.  We  have: 


""  K*  I (n-l)(n+2)(-|J*%.(cosi(») 

a 

^ X (n+2)  (n-1)  (-1  J'*P.  (COS  tl>) 

n 

Hence,  for  r = ro  : 


where 

^ (OOS  ♦) 


We  notice  that  the  relative  error  of  S Ag  is  of  the  same  order  of  magnitude 

3r 

as  that  for  T.  Tims  the  spherical  harmonic  series  for  9Ag  is  convergent  for 

9 r 

R = ro  (cf.  the  gravity  anomalies).  For  this  example,  we  have  the  relation: 


(6gr)n  = 


+ 6T 

a “ 


1 


4.  The  Error  of  the  Potential 


We  are  now  going  to  estimate  the  errors  of  extending  the  potential 
series  (2. 1)  to  the  surface  of  the  earth.  Some  error  estimates  of  this  type 
were  given  by  Cook  (1967)  and  Levallois  (1969).  See  formulae  (4.5  a-b)  and 
below.  See  also  section  4.3 . 
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Minimum  sphere 


Figure  4. 1.  The  potential  is  computed  at  a point  P of  distance  r 
from  the  earth's  center.  The  sphere  of  radius  r is 
denoted  the  sphere  of  computation.  The  topographical 
masses  between  this  sphere  and  the  minimum  sphere 
makes  the  spherical  harmonic  e:q[>ansion  Va  divergent 
at  P. 
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Tbe  Newtonian  potential  of  the  earth  at  an  arbitrary  point  is: 


(4.1) 

where 


d u 


volume  of  the  earth 
Gp 

Newton's  constant  of  gravitation 

density  of  mass 

(ri®  + r®  - 2ri  r cos  0 

radius  of  the  current  point  inside  the  earth 
radius  of  the  point  of  computation  (See  Figure  4. 1) 
geocentric  ang^e  behveen  the  vectors  Fi  and  r 


At  points  outside  the  "minimum  sphere"  (r  > R)  formula  (4. 1)  can  be  expanded 
into  the  following  series: 


(4.2)  v=  JJJ  P„  (COS0)  dl/ 

n = 0 V 


The  corresponding  convergent  series  for  points  inside  the  minimum  sphere 
(r  < R)  is: 


(4.3) 


V,  = 


where  r,  is  the  radius  of  a current  point  at  the  surfoce  of  the  earth 
(r  ^ ri  ^ rt ) and  O is  the  unit  Sfdiere.  See  Figure  4. 1. 

I The  error  fiv  of  extending  formula  (4.2)  to  a point  inside  the  minimum 

sphere  is  given  by  the  difference  between  (4.2)  and  (4.3): 

f (4.4)  6v  6v„ 

n = 0 


i 

I 

f 


where 


(4.5a)  6v„  = 1 JJ  J M [(f  J-  P„  (cos  V*)  d »/  , r.  sr 

a ' 


(4.5b)  «V.  = i JJ  /„  (0“  [(aj”'.  p.  ,o„s  *)  dP 

a 

Formiilae  (4. 5a  -b)  were  derived  by  Cook  (1967)  and  Levallois  (1969).  Cook 
drew  die  conclusion  that  this  error  is  in  the  order  of  Jg  where  Jg  is  the  n-th 
zonal  harmonic  of  the  earth's  gravity  field.  This  result  can  easily  be  combined 
with  Kaula's  rule  of  thumb  (1.1)  for  estimating  the  error  of  a series  of 
spherical  harmonics.  Even  if  the  series  were  extended  to  infinity,  the  error 
would  be  ne^igible,  so  that  " the  'satellite  geoid'  is  a close  enou^  approximation 
to  the  true  geoid".  However,  it  is  not  at  all  obvious  that  6Vg  according  to 
(4.5  a-b)  is  of  order  J®.  Levallois  (1969)  came  to  a different  result,  which 
is  reported  in  section  4.3. 

In  this  section  and  section  4. 1 we  are  going  to  develop  a formula  for 
6 Vg  which  can  be  used  in  a numerical  integration.  From  Heiskanen  and 
Moritz,  1967,  p.  33  we  obtain: 

n 

(4.6)  Pg (cos$)  = y [Rn.(e,X)Rg.(ej,x,)+  Sg.(e,X)Sg.(e,,x, )] 

■=o 

where 

Rn» » Sg,  = fully  normalized  spherical  harmonics. 

We  have 

(4.6a)  JJ  Rg,  SgV  da  = 0 for  all  n,  n"^  m and  m^ 

and 

(4.6b)  1/4  TT  JJ  Ug.  UgV  da  = 6gg' 

where  Ug,  is  any  of  the  harmonics  Rg,  or  Sg,  and  6 is  the  Kronecker's  delta. 
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From  formulae  (4. Sa)  and  (4. 6a  - b)  we  obtain: 


<4.7«) 

where 

(4.7b) 


a 

6V„  = ^ [ a„.  R„.  (8,X)  +’B„  (0,X)] 


■ sO 


V — 1 — r r u r,-  fi . ‘ } 'S>' 


If  we  asaume  that  n = n (8»X)  (independent  of  r),  we  obtain: 


,a, 


where 


(4.8b)  I( 


0 if  r Ti 


{ 


n + 3 


n - 2 


(r,/r)“-l 


- foi  (r,/r) 


if  r < ri , n = 2 


Formulae  (4. 8a  -b)  can  be  applied  with  ti=  Ho=  constant  if  we  neglect  the 
ellipticity  of  die  earth  (spherical  earth  with  topograidiy).  The  integration  is 
then  performed  over  the  topographical  masses.  If  we  also  consider  the  earth's 
ellipticity  the  integration  will  include  the  masses  of  die  oceans  above  the  lower 
bound  r of  the  integration.  See  Figure  4.2.  If  we  assume  that  die  oceans  have 
the  constant  density  ft,  and  the  solid  crust  the  constant  density  Ho  the  int^ral 
(4.7b)  becomes  (see  Figure  4.3): 


(4.8c) 


I pp  — ■■  (•Rb«(®»X)'> 

{b„}"  2^11 " 


where 


Ml(r,r,)  = MoI(r,ii,)  + I(rb,r,) 
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Jd 


1 


Minimum  sphere 


Earth  ellipsoid 
Topography 


Ocean 


Figure  4.2.  The  computation  of  the  error  coefficients  a„  and  bna 

for  any  point  P on  a selected  sphere  of  radius  r (sphere 
of  computation)  includes  the  integration  of  all  masses  of 
the  oceans  and  topography  between  the  spheres  of  radii 
r and  R. 


Figure  4. 3.  The  density  of  the  masses  above  the  sphere  of  computation 
is  iio  for  the  solid  crust  and  Uh  for  the  ocean. 


i(x.y) 


rb  = geocentric  radius  of  tbe  bottom  of  the  ocean. 


The  integral  I (x^y)  is  straight  forward  and  the  result  is  (cf.  4.8b): 


{ 


0 if  X S y 


l(x.y)  = I (y/rr-jx/r)»f.  , i.t6rf.~I-^l.r/x£:\ 


ll/rf-{x/rf  . On  iy/x)  , X <y  , n = 2 
5 


The  coefficients  aj,,  and  can  be  used  to  correct  the  spherical  harmonic 
development  of  the  potential'^at  all  points  on  the  sphere  of  computation  with 
radius  r.  Then  we  add  a correction: 


AC. 


n + l 


to  each  potential  spherical  harmonic  coefficient  of  the  series  exi>ansion.  It 
should  be  emphasized  that  these  corrections  are  valid  only  for  the  potential 
and  cannot  be  used  for  improving  near  surface  expansions  of  gravity  anomalies 
(cf.  section  6). 

From  (4.5  a)  we  draw  the  conclusion  that  if  the  masses  between  the 
spheres  of  radii  r and  R were  symmetrically  distributed  with  respect  to  the 
earth’s  axis  of  rotation,  then  all  tesseral  harmonics  a,,  and  bn.  (with  m^O) 
would  vanish.  This  is  obvious  because  in  that  case  (4. 7b)becomes: 


TT  an 

(4.9) 

S’!!*  ~ r Q n r Rni 

dX  sin  6 

d e 

e=o  x=o 

where 

n+l  J 

(2n+l)r  ^ r 

m+i- 

1 J dr, 

* See  formula  (2.1). 


I: 
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and 


an 

J R„,  dX  = 0 for  m 0 
X=0 


The  same  is  true  for  bn,  . However,  for  the  real  earth  with  its  irregular  mass 
distributions  this  assumption  is  not  valid. 


4.1  An  Approximate  Formula 


Formula  (4. 8b)  can  be  approximated  in  the  following  way  (for 
H = r,  - r > 0): 


k=  0 


1 + /n+3)  ^ fiLt+  (P+3)(n-f2)(n+l)  /H  f , (n^-3)(n+2)(n+l)n  /Hf  , 

' 'r  2 \t  J 1x2x3  \r/  2x3x4  \r  / 


and 


k = 0 

= 1 - (n-2)  ii  + fHt-  jn^g^n-ljn  / H\%  (n-2)(n-l)nm-fl) 

' ' r 2 \r  / 1x2x3  \r/  2x3x4  \r  J 

Inserting  these  series  expansions  into  (4.8b)  we  obtain: 

,4.10,  I,r.r.,  = ,2.«,  H»[i4i 


In  Table  4.1  we  compare  formula  (4.8b)  with  the  approximation 
(2n  + 1)  H®/2  according  to  formula  (4. 10).  It  is  shown  that  the  error  is  increasing 
with  the  degree  (n).  For  n = 50  the  error  is  less  than  1%  and  for  n = 100  it  is 
almost  2%. 
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Table  4. 1 


Comparison  between  (2n+l)  H®/2  and  I(r,R) 
R = 6384  km,  H = R-r=  27  km.  Units  : km^ 


n 

(2n+l)  H®/2 

I(r,R) 

Diff. 

Diff. /I  % 

10 

7654. 5 

7677.4 

22.9 

0.3 

30 

22234. 5 

22328.6 

94.1 

0.4 

50 

36814. 5 

37060.1 

245.6 

0.7 

100 

73264. 5 

74591 

1327 

1.8 

150 

109714. 5 

113812 

4097 

3.6 

200 

146164. 5 

155627 

9462 

6.1 

250 

182614 

201019 

18404 

9.2 

Inserting  (4. 10)  into  (4. 8a)  and  neglecting  terms  of  order  — and  higher 

r 


we  finally  arrive  at: 


(4.11) 


where 


,0  if  r s r, 

H = \ 

<rt-T  r < r* 


If  we  also  consider  the  ellipticity  of  the  earth  we  obtain  from  (4. 8c): 

{4.11a)  ^ JJ  [(/“o-M«)  H?+  H^] 

where 


fio  = density  of  the  solid  crust 
= density  of  Uie  oceans 
^ H if  continent 
Hb  =\  0 if  ocean  with  r ^ rb 
rb-r  otherwise. 
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4.2  Effect  of  Ellipticity 


In  this  section  we  are  going  to  study  the  error  of  extending  the  spherical 
harmonic  expansion  of  the  potential  to  a point  P at  the  surface  of  an  oblate 
homogeneous  ellipsoid  (see  Figure  4.4).  The  error  is  caused  by  the  masses 
outside  the  sphere  of  computation  of  radius  r. 


Figure  4. 4.  The  downward  continuation  error  at  P is 
obtained  by  integration  over  the  (shaded) 
masses  between  the  spheres  of  radii 
r and  R. 


The  error  coefficients  due  to  the  improper  downward  continuation  can  be 
determined  by  formula  (4.11).  With  notations  according  to  Figure  4.4,  we 
obtain: 


r£ 


1-e  sin^  6 


r = R 


1-e  sin  6e 
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and 


I 


_ ^ 

H = te  - r = (cos^  do  - cos®  0)  + O (e*) 

A 


It  has  already  been  shown  in  (4.9)  that  the  coefficients  1,.  and  'Em  are  0 for 
m ^ 0.  For  m = 0 we  obtain  from  (4. 11)  and  the  above  expression  for  H, 
when  neglecting  terms  of  order  higher  tl^  e^ : 

TT-  00 

(4.12a)  I„o  = ce*  J (cos^  0o  - cos®  0)®  (cos  0)  sin  0 d 0 


where 


c = G p n B®/4 


In  the  special  case  0o=  0 we  have: 


(4.12b)  a,o  = ce*  J (1-  cos®  0)®  Pn  (cos  ®)  sinB 


The  factor  (1  - cos®  0)®  is  readily  rewritten: 


(1  - cos®  0)®  = Co  + C2  Pg  (cos  0)  + C4  P4  (cos  0) 


where 


Co-  — , Cg-- 


21/T  ’ 105 


Inserting  this  expression  into  (4.12b)  and  using  the  orthogonality  property  of  the 
Legendre's  polynomials: 


1 

J P,  (t)  P.  (t)  dt  = { 


2 n = m 
0 n ^ m 


1 


f 


we  finally  arrive  at 


a 0 Q 2 c e 


‘ { 


Co 

if  n = 0 

Ca 

if  n = 2 

C4 

if  n = 4 

0 

otherwise 

For  p = 2.67g/cm®,  R = 6378  km,  ro  = 6371  km,  y = 978  gal  and 


e = 0.0067  we  obtain: 


27.9  m 

y 


^=-17.8  m 

y 


4.0m 


^1  4.3  X 10'® 

yro' 


—2^  1 <=«  5. 8 X 10 
yroCgo 


-a 


S' 40 


y ro  C40 


1.1 


From  these  coefficients  we  obtain  the  following  downward  continuation 
error  at  the  poles  ( 0o  = 0): 


6v  = y a„o  '^2n+l  = 


■ = 0 


Thus  the  errors  of  the  individual  spherical  harmonics  compensate  each  other  at 
the  poles. 

For  an  arbitrary  point  P ( So)  on  the  surface  of  the  ellipsoid  the  downward 
continuation  error  in  an  e}q)ansion  to  degree  N is  given  by: 


6 V(N)  = y a„o  (r)  P„  (cos  Oq) 


n s 0 


Inserting  (4. 12a)  into  this  formula  and  changing  the  order  of  summation  and  integration 
we  obtain: 
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Hfi  V * 


ft  V (N)  = c e*  J (COS®  00  - cos®  0f V (cos  0)  P,  (cos  0o)  sin  0 d 0 

00  »=  0 


Using  the  substitutions: 


t = cos  0 , to=  cos  00 


Pn  (t)  Pn  (to)  = (2n+l)  P„  (t)  P„  (to) 


we  arrive  at 


(4.13)  6V(N)  = ce^  J (to®  - t®)®  Y (2n+l)  P„  (t)  P„  (to)  dt 

-to  "=  0 

It  is  shown  in  Proposition  4 of  the  Appendix  that  this  formula  may  be  written: 


6V(N)  =c  e*tb®  gM(f) 


where 


2M->-3 


IV-i  (to)  1 


N if  N is  odd 


/■  U ID  « 

^ N+1  if  N is 


0 s rs  to 


The  function  gM  (and  subsequently  6V(N))  has  the  following  properties: 


Kn  (0)  = gN  (1)  = 0 
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and 


lim  gN  = 0 

N-*oo 


From  formula  (6.9)  we  obtain: 

ce®  gN(to)  = 6^(N) 

and  the  loug^  aiq>roximation: 

(4.13a)  6v  (N)  ^ t|  6 Ag  (N) 

4 

where  5Ag(N)  is  the  downward  continuation  error  of  Ag  in  an  ejqjansion  to  degree  N. 
For  R = 6378  km,  e®  = 0.0067  and  y = 978  gal  we  have: 


0.011  io  6 ^(N)  meter, 

where  6 Ag  is  in  units  of  mgal.  From  Fig.  6.1  we  finally  obtain: 

j 

0.45  m and  0.19  m i 

y y 

This  development  indicates  that  the  downward  continuation  error  of  the  truncated 
spherical  harmonic  expansion  of  V is  small  and  decreasing  with  the  degree  of 
truncation.  On  the  other  hand,  the  following  spherical  harmonic  expansion  for 
the  potential  in  the  exterior  of  a homogeneous  oblate  ell4>soid  is  known  from 
potential  theory  (see  MacMillan,  1958,  p.  363): 


where  M is  the  mass  of  the  ellipsoid.  This  devel(^ment  has  the  radius  of 
convergence  ae,  which  is  inside  the  ellipsoid.  However,  the  potential  provided 
by  (4. 14)  is  correct  only  on  the  surlEace  and  outside  the  ellipsoid.  This  fact  is 
indirectly  verified  by  letting  e approach  zero,  in  which  case  (4. 14)  becomes: 


This  is  the  well  known  exterior  potential  of  a homogeneous  sphere.  The  radius 
of  convergence  is  zero. 

In  conclusion,  formula  (4. 13)  is  the  difference  between  the  expansion 
(4.14)  to  degree  N and  a spherical  harmonic  e}q>ansion  (at  the  radius  r),  which 
is  valid  also  inside  the  surface  of  matter.  For  N approaching  infinity,  the  sums 
of  the  two  series  are  identical  for  exterior  points.  Inside  the  surfkce  of  matter 
they  differ. 


In  geodesy,  we  are  mainly  interested  in  the  exterior  gravity  field  and 
a potential  expansion  of  the  mean  earth  ellipsoid  (MEE)  similar  to  (4.14)  is  there- 
fore most  convenient.  Unfortunately,  this  ellipsoid  is  not  homogeneous  and  (4. 14) 
can  not  be  applied.  However,  the  MEE  is  currently  approximated  by  a level 
ellipsoid,  the  external  gravitational  potential  of  which  can  be  expanded  in  the 
following  way  [Heiskanen  and  Moritz,  1967,  formulae  (2-88)  and  (2-91)]: 


< 


-t-D” 

(2n+])(2n+3) 


n+5n  ) f-^)  Pan(coBd) 

Ma^e=  ^ ^ 


where  C and  A are  the  principal  moments  of  inertia.  If  we  anticipate  this  model 
for  the  mean  earth  ellipsoid,  there  will  be  no  contribution  to  the  downward 
continuation  error  from  the  ell^soid.  However,  the  previously  derived  correction 
formulae  (4. 8 a-b)  and  (4. 11)  for  the  distuibing  topography  needs  a modification. 

In  these  formulae  it  was  assumed  that  the  radius  (r  .)  of  the  mean  sea  level  was 
constant.  Now  we  assume  that: 


r.  = R 


1-e" 

l-e"*  sin"*  e 


Then  (4. 8 a-b)  and  (4. 11)  are  modified  by  inserting: 


(4.15) 


= f r.  if  r<r, 
I r if  r 5 r , 


Formulae  (4.8a-b),  (4.11)  and  (4.15)  should  be  feasible  expressions  for  a 
numerical  study  of  the  downward  continuation  error  of  the  potential.  Formulae 
(4.8c)  and  (4. 11a)  are  less  convenient  unless  the  spherical  harmonic  expansion 
that  is  valid  at  the  entire  sphere  of  computation  (also  inside  the  earth)  is  required. 


4.3  An  Error  Estimate  for  the  Potential 


We  may  substitute  6 V of  formulae  (4.4)  through  (4.7)  by  the  error  of 
the  disturbing  potential,  6t,  if  we  assume  that  the  reference  field  is  unchanged. 
Formula  (4.5b)  is  then  written: 

(4.16)  6T.=  i JJJ  „(0"[(aj’*'.i]  p.,coM)dt- 

a r 


At  sea  level  (r  = ro  = R-H)  this  error  is  of  the  following  order  of  magnitude: 


(4. 17a) 
where 
(4. 17b) 


6T„ 


-[(tr-i 


ATn 


° T,  II  I " *> 

to 


For  low  degrees  of  n we  may  use  the  following  approximation: 

(4.18)  =%*  (2n+l)  AT, 

R 

AT,  is  the  contribution  to  T„  from  the  topography.  If  we  assume  that  the 
disturbing  potential  is  generated  entirely  within  the  topognq>hy  of  the  earth,  we 
have  [cf.  formula  (4.3)]: 

(4.19)  AT,  T„ 


Assuming  that  the  low  order  harmonics  of  T is  only  partly  originated  within  the 
topography,  we  may  use  the  following  approximation: 
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(4.20) 


At,  = PnT, 


where 


P»  = 1 - e' 


0 < a < ® 


The  coefficient  a can  be  determined  empirically  by  conq>aring  a spherical 
harmonic  expansion  of  AT  (formula  4. 17b)  with  that  of  T.  In  the  limit  a=  °° , 
(4.20)  is  identical  with  (4.19).  For  a = 0,  we  obtain  AT,  = 0,  that  is,  no 
contribution  to  T,  from  the  topography.  The  approximations  (4. 18)  and  (4. 19) 
were  used  by  Levallois  (1969)  in  the  following  way: 


6t  = y 6t„  § y (2n+l)  T„ 

a « 


= H y 2H  Ag 


where  we  have  used  the  relation: 


T = ^ Agn 

" n-1 


The  error  of  the  geoidal  height  is  then: 


6n  = — « 2H 


For  H = 27  km,  ^ = 100  mgal  and  y=  980  gal  this  formula  implies: 


6N'=«  5m 


However,  the  ai^roximation  of  (4. 17a)  by  (4. 18)  is  adequate  only  for  lower  orders, 
(say,  n < 150).  Consequently,  the  above  order  of  the  error  does  not  hold  if  we 
include  harmonics  to  infinity. 


i 


4.4  Global  RMS  Errors 


Spherical  harmonics  of  different  degrees  and/or  orders  are  orthogonal 
to  each  other,  see  formulae  (4.6a-b). 


\ • 

( 

! 

Thus,  the  global  RMS  error  of  fiT  is  given  by: 

1 

1 

1 

i 

(4.21a) 

n a 

i 

where 

; 

(4.21b) 

^ If  (®Tn)^  da  = ||6tJP  (C+ 

■ = 0 

j 

1 

The  last  formula  is  obtained  by  substituting  (fiV)n  by  (fiT)n  in  formula  (4.7a). 

The  total  error  of  T for  a truncated  series  of  spherical  harmonics  is  (cf.  section  3) 

(4.22) 

€t  (N)  = 6T(N)  + er  (N) 

where 
ei  (N) 

N is  the  order  of  truncation,  6t  (N)  is  the  downward  continuation  error  and 
is  the  truncation  error.  We  obtain: 

(4.23) 

N+l  N+1 

where 


do 


o^  are  the  degree  variances  of  the  potential.  From  (4.22)  we  arrive  at  the  following  I 

formula  for  the  total  RMS  error  when  noting  the  orthogonality  between  spherical  j 

harmonics  of  different  degrees:  ; 

i 

1 

||ct(N)11®  = I16T(N)||^  lleT(N)||=  ; 


Definition  4.1:  The  optimum  degp^ee  of  truncation  (Nopt)  implies 


II  fT(N)||  = minimum 
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The  downward  continuation  error  is  increasing  with  K wfaUe  the  truncation  error 
is  decreasing.  Consequently,  It  €t  (N)||  is  minimum  when  the  Nth  term  of 
116t(N)H  equals  the  (N+l)th  term  of  Her  (N)  11,  i.e. ; 


An  approximation  of  6Tn  is  obtained  from  (4. 17a)  and  (4. 19) 


where  R is  the  radius  of  the  minimum  sphere  (~  6384  km)  and  ro  is  the  mean  earth 
radius  (~6370  km).  Inserting  (4.25)  into  (4.24)  we  obtain: 


.aN+i 


As  the  degree  variances  are  decreasing  with  N we  easily  arrive  at; 


or  with  the  above  numerical  values  for  R and  ro 


The  following  values  for  the  potential  degree  variances  were  given  by  Tscheming 
and  Rapp  (1974): 


(n-2)(n-l)(n+24) 


where 


(TBja/ro)  = 0.999617 
425.18  [rngal®! 

6371  [km]  (ro  = 6370  km  was  used  in  the  computations.) 


Figure  4. 5.  The  RMS  errors  of  the  downward  continuation  and  truncation  of  the  geuidal 
undulation  expansion  to  degree  N.  R = 6384  km  and  r©  = 6370  km. 


A 


The  graphs  l|6T<N)ll/y,  lleT(N)ll/y  and  llcr(N)l|/y  for  these  degree  variances 
are  illustrated  in  Figure  4.5.  The  total  error  H €t  ||/y  is  minimum  for  N = 157. 
This  error  is  1. 17  m («  4. 7%  relative  error).  Most  of  this  error  (1. 05m)  is  due 
to  the  downward  continuation  error. 

Ui  the  estimates  above  we  have  assumed  that: 

r, 

T„  ^ AT„  = i J|  J H Pn  (cos  >l>)  dv 

a To 

A different  aiq>roximation  was  suggested  in  formula  (4.20): 


AT„  = p„T„ 

where 

p„=  l-e“®’  , 0<  a < » 


I 


I 


I 


For  this  estimate  (4.25)  and  (4.26)  take  the  forms: 

(4.25')  6t„=[(^J  -l]pnT„ 

and 

ro^(l+^) 

(4.26)  Nopt  ^ 


Let  us  assume  that : 


Pioo  = 0.95 

i-ioo 

This  relaticm  corresponds  to  a = 0.03.  The  resulting  RMS  errors  for  downward 
continuaUonare  illustrated  in  Figure  4.5.  It  is  shown  that  ||6t||  < 4%  of  ||t||  for 
N < 200  and  the  total  RMS  error  is  minimum  for  N = 158.  This  minimum  is 
approximately  3.5%  of  11t|1.  The  largest  uncertainty  of  these  estimates  is  the 
relation  _R_  introduced  in  the  approximation  (4.17a).  The  sensitivity  of  the 
To 

solutions  of  116t11  and  Nopt  for  the  difference  H = R-ro  is  obtained  from  the 
following  approximate  solutions: 


1 


I 
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6T(N)  -S  Y (2n+l)  T„ 
R 

n!=0 


ro^2 


Differentiating  this  formulae  with  respect  to  H we  arrive  at: 


(4.27) 


|d6T(N)|_  idNopt  I = 1^1 
' 6T(N)'  ' Nopt  H ' 


Now  H is  in  the  order  of  13  km  and  the  ^ror  dH  may  be  several  kilometers. 

If  we  assume  that  dH  4 km,  we  obtain  a relative  error  of  31%.  Thus  we  conclude 
that  the  RMS  estimates  given  in  this  section  mi^t  have  considerable  errors  due 
to  the  uncertainty  of  the  basic  approximation  (4. 17a). 


5.  The  Error  of  Gravity  Disturbances 


The  error  of  downward  continuation  of  the  vertical  derivative  of  V in 
a series  expansion  is  obtained  from  formulae  (4.2)  through  (4.4): 


6 ^ ^ ^ = -L  6V 

dr  dr  dr  dr 


If  we  subtract  the  normal  field  U from  V in  this  differentiation  (and  change  sign 
of  the  whole  e^qpression)  we  obtain  the  gravity  disturbance: 


Assuming  that  the  normal  field  is  the  same  for  T#(=V,-U)  and  Ti  (=Vi-  U)  we 
have  the  error: 


6t  =-^6v 

dr  dr 


where  according  to  formula  (4. 4)  and  (4. 5a): 


“ T Ji  I " [(f*)'-  (f  ri  ** 


1=0  O r 


We  use  the  notation: 


A6r  =7  4, 


wheie 


4„  = 


n=0 


It  is  shown  in  the  Appendix^  Proposition  1,  that  An  can  be  written: 

JI J " [(^J  * iB  "] 


o r 


This  expression  can  be  rearranged  in  the  following  way: 

r. 


A„  = 


n+1 

^n+3 


2n+l. 


a r 


Comparing  with  formula  (4. 5a)  we  find: 


(5.2)  A„=  6t„  + ^ JJ  J M (0  P„(cosV))  dP 
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where  6 T„  is  the  error  of  T„  caused  by  the  analytic  continuation  down  to  r. 
Furthermore,  we  have  for  r = ro  [radius  at  mean  sea  level,  cf.  (4. 17b)]: 


^ If  J " {fj 


To  M J \ ri , 

a To 


AT„  is  the  contribution  to  the  disturbing  potential  harmonic  Tn  from  the  anomalous 
masses  located  outside  the  sphere  of  radius  ro.  Thus  we  have: 


A„=  — [(n+1)  6t„  + (2n+l)  AT„] 
To 


or,  after  inserting  (4. 17a)  for  6Tn: 


Finally,  we  sum  the  harmonics  A „ to  the  degree  of  truncation  (N): 


where  according  to  (4.20): 


^Tn  = Pj  Tn  , Pn  — 1 ~ e 


Formula  (5.2)  gives  an  "exact"  value  of  the  downward  continuation  error  for 
3 T 

- _ Formula  (5. 5)  is  an  approximation. 


6.  The  Error  of  Gravity  Anomalies 


The  gravity  anomalies  are  related  to  the  gravity  disturbances  by: 
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Ag=  6, 


2T 

r 


III  the  same  way  we  obtain  for  the  errors  of  downward  continuation: 


6Ag  = A6,  - 

r 

or 

6Ag„  = A„- 

r 


where  is  the  error  of  and  is  its  n-th  harmonic.  By  inserting 

(5.2)  we  obtain: 

(6.1)  6Ag.  = ^ 6t„  + JJ  J ^ 

o r ‘ 


and  approximately  from  (5.4)  for  r = r©  : 


(6.1')  6Ag„=  [^(n-1)  +n  + 2] 

This  error  is  of  the  same  order  as  An  . 

The  first  part  of  (6.1)  is  already  expressed  in  terms  of  spherical 
harmonics  (see  formulae  4.7  a-b).  The  second  term  may  be  written; 

r< 


(6.2a) 

n 

— 7 IC„.  R..+  Jn,  S,.] 

T*  _i 

^ a r 

* A 

■ sO 

where 

if:)  - 

(6.2b) 

5: 

O=UJ<‘(0  • 

Assuming  that  n = tnd,X)  (independent  of  r)  we  obtain: 

i 

I 

j 


1 

I 


) 


t 


I 
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(6.3a) 


where 


(6.3b) 


(aj'  If" 


J(r.r,)=r®( 

C n-2 

^(n/r) 


,B-a 


, If  r a r, 

, If  r < r,  , n ^ 2 
, If  r < r,  , n = 2 


For  H = r.-r>0  the  last  formula  can  be  expanded  into  the  following  series: 


(6.3c) 


J(r,  r.)  = rH  - if  + ... 


The  error  of  ^ is  thus  given  by: 

00  0 

(6.4)  6/^  = y y [A„  R„  + Bn,  S„,] 

ttssO  BsO 

where 


or 

(6.5) 

where 


^n,  [(0~1)  ®n,  Cn,  ] 


Bn,  ^ [(**~1)  "^n,  + dn,  ] 


- ■ 


if  r 2 r. 


K(r,r.)  = r ■[  (n-1)  ^-*"^{3  - (n+2)  if  r<r„n?f2 


+ 4 (r,/r)  if  r < r,  , n=  2 
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Tbe  last  formula  may  be  approximated  by: 


(6.6) 


lb 


For  one  more  term  in  the  approximation  (6. 6)  see  Prop.  2 of  the  Appendix. 
Formula  (6.6)  can  also  be  determined  directly  from  the  potential  coefficients 
Sna  and  baa*  using  the  boundary  condition  (2.2): 


or 


K(r.r.)  = -(|^  + j)  I(r.r.) 


Inserting  (4. 10)  and  noting  that  ~ ~ we  finally  obtain: 

or  oH 


(6.7) 


K (r.r.)  = (2n+l)H  [l  + + •••  ] 


where  the  first  term  of  K is  the  same  as  in  (6.6).  In  Table  6. 1 it  is  shown  that 
(6.6)  is  a good  approximation  of  (6.5). 


Table  6. 1 


Comparison  between  H = 27  km  and  - — K (r,  R) 

2n+l 


for  R = 6384  km,  r=R-H 


10 

20 

50 

100 


2II+1 


K(r,  R)  [km] 


27.00877 

27.03394 

27.20732 

27.82721 


Difference  [m] 


- 8.77 

- 33.94 
-207. 3 
-827. 2 


Rel.  diff.  % 


0.03 

0.13 

0.77 

3.06 


-38- 


In  formulae  (6. 3 a-b)  and  (6.6)  we  have  assumed  a spherical  earth  with 
topography,  where  the  only  contribution  to  the  coefficients  An.  and  Bu  are  given 
by  the  topography.  We  may  also  include  the  effect  of  the  earth's  ellipticity.  In 
that  case  ^K(r,  r.)  of  (6.5)  is  replaced  by  (see  Figure  4.3): 


Ho  K(r,rt,)  + K (rt,  r,) 


where  rb  is  the  radius  of  the  ocean  bottom  and: 


0 if  rb  2s  r. 


:(rb,r,)=  r j (n-1) 


’ n#  2 


f-(r^/rf  _ , i\  < r, 

5 +4«,(r,/r.)  . 


In  this  case  we  arrive  at  the  following  approximate  formula  corresponding  to  (6.6); 


(6.6  ) 


where 


F(H,Hb)  = M-  H+  (Mo-M«)  Hb 


H if  continent 

Hb  = I 0 if  ocean  with  r > rb 


rb  - r otherwise 


However,  in  the  next  section  it  is  shown  diat  formula  (6. 6 ) is  not  very  usefel  for 
our  purpos'  k>  determine  the  downward  extension  error  to  the  surface  of  the  earth. 
It  seems  more  convenient  to  use  a modification  of  (6.6)  [cf.  section  4.2]. 
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6. 1 Effect  of  Elliptic  it 


The  error  of  Ag  due  to  the  improper  downward  continuation  to  the 
surface  of  a homogeneous  ellipsoid  can  be  estimated  in  the  same  way  as  the 
error  of  V in  section  4.2.  In  this  case,  we  obtain  for  an  arbitrary  point 
P(6o)  on  the  ellipsoid. 


where 


6 Ag  (N)  = y A„o  Pn  (cos  00) 


Sbo  = ke^  J (008^00  - cos^  0)  P„  (cos  0)  sin  0 d 0 


k = GPo  tt  R 


For  00  = 0 the  last  formula  is  easily  integrated  (cf.  section  4.2)  and  the  result  is; 


-y  for  n = 0 


A„o  = 4ke* 


I -xJl 


0 otherwise 


For  Po  = 2.67  g/cm®,  R = 637S  km  and  e^  = 0.0067  we  have  the  following  numerical 
values  of  the  error  coefficients: 


Aoo  — 3. 188  gal. 


Sao  = - 1.425  gal. 


Inserting  these  coefficients  into  (6.8)  we  arrive  at: 


6Ag(N) = y Abo  Pb  (1)  = 0,  N i 2 


4 


I 


i 

I 

I 

I 


Thus  there  is  no  error  at  the  poles.  For  an  arbitrary  point  P(0o)  we  obtain  from 
Proposition  3 of  the  Appendix: 


(6.9) 

where 


6Ag(N)=  2ke=  ft-.  ] 

Pii  = Rn  (cos  % 

, N if  N is  odd 

M = i 

l-N+1  if  N is  even 


From  (6.9)  it  follows  that  6Ag  (N)  is  decreasing  towards  zero  for  increasing  N, 
because: 


R 0 when  N - ® 


For  low  order  e;)q}ansions  it  is  shown  in  Figure  6. 1 that  the  "errors” 
according  to  (6.9)  are  considerable.  For  example,  the  maximum  error  for 
N = 16  is  83  mgal.  Again,  it  should  be  emphasized  that  these  values  are  the 
deviations  of  the  external  series  from  the  true  series  at  the  sphere  of  computa- 
tion (cf.  section  4.2).  When  N approaches  infinity  the  sums  of  these  two  series 
are  identical  for  exterior  points.  Thus  the  external  series  is  indeed  convergent 
and  is  a correct  representation  all  the  way  down  to  the  surface  of  a homogeneous 
ellipsoid.  The  same  holds  for  a level  ellipsoid.  Hence,  if  we  anticipate  the 
mean  earth  ellipsoid  as  a level  ellipsoid,  the  only  contribution  to  the  downward 
continuation  error  of  Ag  to  the  earth  is  given  by  the  topography.  The  error 
coefficients  A„,  and  Bn«  are  then  computed  from  formula  (6.5)  or  (6.6)  with 
the  modification  (cf.  section  4.2): 

(-r.  if  r<  r. 

r = 1 

*«  r if  r s r, 

r,  = Ry(l-e®)/(l-e^sin®0) 


Formula  (6.6  ) should  be  used  if  the  spherical  harmonic  expansion  that  is  valid 
or  Ihe  entire  sphere  of  computation  is  required. 


1 


1 
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The  downward  continuation  "error"  of  a 


r 

f 

i 

[ 

\ 

f 

i 

\ 


i 

I 

i 


I 

1 

f- 


6. 2 Global  RMS  Errors 


The  total  error  of  for  a truncated  series  of  spherical  harmonics  is: 


(N)  = 6 Ag  (N)  + e^  (N) 


where  6^  (N)  is  the  error  for  the  continuation  and  e^  (N)  is  the  truncation 
error.  The  global  RMS  of  (N)  is: 


(6.10)  1U/).(N)1|=  [ll6Ag(N)f  + lle^(N)||^l^ 


At  sea  level  we  have ; 

00 

lle^.(N)ll^  = y (Ag) 

N+l 


or,  for 


(Ag)  = ^ 
To 


(6.11) 


l|eA.(N)ll^  = y (^) 

N+l 


where  are  the  disturbing  potential  degree  variances. 

From  (6.  l')  and  (4.20)  we  finally  obtain  the  approximate  formula: 


(6.12)  |l6Ag(N)||^  " ^ i 


n=0 


where 


Pn  = 1 - e 


-Q(b 


a = 0.03 
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The  formulae  (6. 10)  - (6. 12)  are  shown  in  Figure  6. 2 for  ihe  degree  variances 
of  Tscheming  and  Rapp  (1974).  In  Ihis  figure  we  also  depict  the  corresponding 
RMS  errors  for  5°x  5°  block  mean  anomalies,  which  are  obtained  by  multiplying 
each  degree  variance  of  formulae  (6. 10)  - (6. 12)  by  the  smoothing  factor  ^n^ 
(see  section  A.  3 of  the  Appendix  or  Meissl,  1971). 


7.  The  Error  of  Vertical  Gradients  of  Gray  it: 


The  error  of  the  vertical  gradient  of  gravity  is  given  by  the  derivative 
of  6 Ag  with  respect  to  r: 


6Ag=y  ± 6Ag, 

n=0 


From  formula  (6.1)  we  obtain: 


(7.2)  A 6Ag„  = - ^ 6T„  + Sfi  ^ 6t„  - AT„  + ^ ^ AT, 


r®  " r 9 r 


where  AT„  is  given  in  formula  (4. 17b).  The  derivative  of  AT„  with  respect  to 
r is  (cf.  Proposition  1,  Appendix): 


^ ATn  = AT„  - r JJ  m Pn  (cos  0)  d o 


Using  the  approximation  yL=  Ho  = constant,  the  last  integral  vanishes  for  n 0. 
Thus  we  have: 


^ AT„  — AT,  , n?^0 
or  r 


Furtheimore,  we  obtain  from  (5.2): 


^ 6T.  = - ^ 6t.  - <2nrt)  ^ 
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Inserting  these  derivatives  into  (7. 2)  we  arrive  at  the  following  formula; 


This  relation  is  the  same  as  the  relation  between  the  spherical  harmonics 
and  Tb  themselves  in  the  external  case.  Thus  we  have; 


The  error  of  the  vertical  gradient  of  gravity  can  also  be  estimated  from  (7. 1) 
(6.4),  (6.5)  and  (6.7).  From  (6.4)  we  obtain: 


where 


Using  the  approximations  (6. 5)  and  (6. 7)  we  arrive  at; 


(2n+l) 


Comparing  this  formula  with  (4. 11)  we  finally  obtain 


° hr 

6T, 

SAg„ 

T„ 

9r 

[■' 


rC..,  _ (n-l)(p-h2) 

^D..J  r 


This  formula  is  equivalent  with  (7. 3). 


8.  Computations  with  5°x  5°  Mean  Elevations 


1654  mean  elevation  blocks  were  used  to  compute  the  downward 
continuation  errors  in  a spherical  harmonic  expansion  of  the  disturbing  potential 
and  the  gravity  anomaly.  In  all  computations,  the  density  of  mass  ( Po ) was  set 
to  2.67  g/cm^.  The  spherical  harmonic  series  were  expanded  to  degree  16  (if 
not  specified). 


8. 1 Computations  for  a Spherical  Mean  Earth 

In  the  first  set  of  computations  we  assumed  a spherical  mean  earth 
with  radius  r = 6371  km.  The  downward  continuation  error  for  each  harmonic 
of  the  disturbing  potential  and  gravity  anomaly  can  then  be  written  according  to 
formulae  (4.11)  and  (6.6)  respectively: 


o 

, r,  - r if  r,  > r 
H = ] 

0 otherwise 

The  computed  error  coefficients  are  shown  in  Tables  A.l  - A. 2.  From  these 
coefficients  the  errors  of  V and  Ag  were  computed  from: 

la  n 

(8.3)  6v=yy  [an.Rn.(e,A)+ b^S„(0,X)] 


(8.1) 

and 

(8.2) 

where 


3 


1 
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L-t 


1 


and 

le  n 

(8-4)  6Ag  = y y [A„,  R„,  (0,A)  + B„  S^.  (6, A)] 

-•J 

nsO  asO 

The  only  area  with  any  significant  errors  6 v/y  was  souUi  of  the  Himalayas.  In 
order  to  estimate  the  errors  at  the  surface  of  the  earth  in  this  area,  a second 
computation  was  performed  with  r = 6374  km.  Then  the  errors  were  interpolated 
between  the  computations  with  r = 6371  km  and  r = 6374  km.  The  maximum  error 
estimate  did  not  exceed  half  a meter. 

The  errors  of  the  gravity  anomalies  were  in  the  same  order  as  the 
anomalies  themselves  (Fig.  8.1).  These  surprisii^ly  large  errors  are  not  in 
agreement  with  the  empirical  results  obtained  through  a direct  comparison  of 
satellite  derived  spherical  harmonic  expansions  with  terrestrial  gravity  anomalies. 
We  conclude  that  the  approximation  of  the  mean  earth  with  a sphere  is  too  rough 
for  estimating  the  anomaly  errors. 


Figure  8. 1.  The  downward  continuation  error  of  Ag  along  two  profiles 
at  sea  level  using  formulae  (8. 2)  and  (8. 4).  Degree  of 
truncation:  16. 
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8.2  Computations  for  an  Ellipsoidal  Mean  Earth 


We  anticipated  an  ellipsoidal  mean  earth  with  the  following  parameters. 


a = 6378140  km 


f = 1/298.257 


The  integration  of  the  error  coefficients  for  6v  and  6^  was  performed  according 
to  formulae  (4.11a)  and  (6.6^): 


(sll  = ° IJ  + H) 


where  Po  = 2.67  g/cm^  and  Pw  = 1.03  g/cm® 


The  computed  coefficients  aoo,  ago  . a4o  , Aoo  . and  Ago  agreed  very 
well  with  those  determined  in  sections  4. 2 and  6. 1.  The  downward  continuation 
"errors"  for  two  profiles  are  shown  in  Figures  8.2  and  8.3.  It  was  found  that  the 
geoidal  undulation  errors  were  small  (less  than  a meter)  while  the  gravity  anomalies 
deviated  up  to  70  mgal  in  an  expansion  to  degree  16  (cf.  sections  4.2  and  6. 1).  The 
results  are  not  in  agreement  with  the  empirical  knowledge  of  the  errors  of  down- 
ward continuation. 


8. 3 Computations  for  a Level  Ellipsoid  with  Topof 


In  this  section  we  make  use  of  a mean  earth  level  ellipsoid  with  the 
same  dimensions  as  the  ellipsoid  of  the  previous  section.  For  a level  ellipsoid 
there  i.c  no  contribution  to  the  downward  continuation  error  to  the  surface.  We 
could  therefore  use  the  formulae  (8.1)  - (8.4)  with  the  modification; 


_r  r,  if  r,  > r 

^ r if  r,  s r 
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Figure  8.6.  Downward  continuatioa  errors 
for  ^ according  to  formulae  , 
(8.2),  (8.4)  and  (8.7).  N > 16. 


<^= -67.5®,  r= 6360  km 
(i=-l2.51r=6377  km 


where 


f(l-e^)/(l-e®  sin^  0) 


For  r = a/l-  e®  (at  the  poles),  these  computations  are  identical  with  those  given 
in  section  8.1.  In  Tables  A.  3 - A.  4,  we  give  the  error  coefficients  for  r = 6371  km. 
In  Figure  8.4,  we  show  the  degree  variances  (af)  of  the  potential  coefficients  as 
defined  by: 


(gm)  t 


In  Figure  8.5  the  accumulated  RMS  errors  of  the  geoidal  undulations  to  degree 
N„x  = 16  are  shown.  These  errors  do  not  exceed  0.13  m. 

The  downward  continuation  errors  of  the  gravity  anomalies  in  an 
expansion  to  degree  16  are  depicted  in  Figure  8. 6 along  two  latitudes  at  sea  level. 
In  Figure  8.7,  the  expansions  to  N = 16  and  N = 24  are  compared  for  a profile. 
Large  error  estimates  are  obtained  at  the  edges  of  the  continents.  Otherwise, 
the  errors  are  within  ± 5 mgal. 


9.  Conclusions  and  Final  Remaiks 


In  this  report  we  have  investigated  the  downward  continuation  errors 
for  spherical  harmonic  expansions  of  the  gravity  field  of  the  earth.  The  simple 
model  in  section  3 showed  that  these  errors  are  increasing  with  the  height  of  the 
surrounding  topography  above  the  computation  point.  Furthermore,  the  relative 
errors  are  increasing  with  the  degree  of  truncation  of  the  series.  It  is  also 
shown  that  the  usual  relation  between  the  spherical  harmonics  of  gravity  anomalies 
and  disturbing  potentials  is  not  valid  for  the  errors  of  downward  continuation: 


6Ag„9^  Sll  6t. 


The  errors  of  Agn  are  usually  more  serious  than  would  be  the  case  if 
this  relation  were  true.  One  reason  for  this  is  that  while  the  potential  is  dependent 
only  upon  the  distance  to  the  generating  masses,  the  influence  on  the  gravity 
anomalies  from  masses  located  above  the  computation  point  usually  has  opposite 
sign  to  that  obtained  in  the  downward  continuation  procedure.  Thus,  the  relative 
errors  of  the  gravity  spherical  harmonics  are  more  than  100%  in  the  model  of 
section  3,  where  the  only  contribution  to  die  anomaly  is  the  disturbing  mass  out- 
side the  main  sphere. 


r 


i 

[ For  the  error  of  the  vertical  gradient  of  gravity  (6gr)f  the  following 

t formula  was  found  for  the  model: 


(9.2) 


(6gr)n 


6T 

3 ^ -I  n 


which  is  the  same  as  the  relation  between  the  harmonics  themselves. 


Some  general  global  RMS  errors  were  estimated  based  on  Tscheming/ 
Rapp's  degree  variances.  In  the  low  order  expansions  (N  = 30)  the  errors  of 
geoidal  undulations  are  2.2  dm  in  the  most  optimistic  estimate  (a  = 0.03)  while 
the  RMS  errors  of  gravity  anomalies  reach  16  mgal.  These  error  estimates  of 
the  gravity  anomalies  seem  too  large,  which  might  be  due  to  the  uncertainty  in 
the  basic  approximation  (4.17a). 

Formulae  were  developed  for  a numerical  integration  of  the  downward 
continuation  errors  [formulae  (8.1)  through  (8.7)].  Formulae  (9.1)  and  (9.2) 
were  found  to  be  valid  also  for  the  real  earth. 


In  sections  4.2  and  6. 1 formulae  (8. 1)  and  (8.2)  were  tested  at  the 
surface  of  a homogeneous  ellipsoid.  The  calculations  showed  that  the  "error" 
coefficients  of  degree  0,  2 and  4 were  considerable.  However,  the  errors 
6 V (N)  and  6^(N)  were  attenuating  with  N to  zero  for  N approaching  infinity. 
The  reason  for  this  strange  result  is  the  following.  Inside  the  minimum  sphere, 
the  potential  of  a homogeneous  ellipsoid  may  be  represented  by  two  different 
spherical  harmonic  series:  one  which  is  convergent  on  the  entire  sphere  of 
computation  (also  inside  the  ellipsoid)  and  one  which  is  valid  only  outside  the 
surface  of  matter.  The  latter  is  identical  with  the  expansion  outside  the  bounding 
sphere.  As  we  are  only  interested  in  the  errors  on  and  outside  the  surface,  we 
are  looking  for  the  downward  continuation  error  of  the  latter  development,  while 
formulae  (8.1)  and  (8.2)  give  the  errors  relative  to  the  first  series. 


It  was  therefore  suggested  that  formulae  (8. 1)  and  (8.2)  should  be 
modified  (when  applied  to  the  real  earth)  in  such  a way  that  there  are  no  con- 
tributions to  the  downward  extension  errors  from  the  masses  of  the  mean  earth 
ellipsoid  (level  ellipsoid).  This  modification  [formula  (8. 7)]  implies  that  the 
errors  of  the  spherical  harmonic  expansions  are  merely  due  to  the  topographical 
masses  above  the  sphere  of  computation. 

In  the  final  computations  with  1654  5°x  5°  mean  elevations  (section  8.3) 
we  arrived  at  very  small  errors  of  the  geoidal  undulations  (RMS  error  = 0. 13  m 
for  N = 16).  The  degpree  variances  of  v were  found  to  have  a maximum  38  x 10"^® 
for  n = 6 and  r = 6356.73  km  (at  the  poles). 


I 
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Finally,  we  conclude  ifaat  the  gravity  anomaly  errors  seem  generally 
to  be  within  ± 5 mi^.  However,  at  the  edges  of  the  continents  and  within  rough 
areas  on  the  continents,  larger  errors  might  be  expected.  Such  areas  may  be 
studied  in  detail  by  utilizing  the  method  described  in  section  8.3. 
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Appendix  A.  2 


A.  2 Derivation  of  Some  Formulae 


Proposition  1: 


(A.l) 


where 


= 6v, 


‘ “ HI  J “ [ - (i)”'  ] '■■  ^ 


implies 


d F = ri^  d rj  d O' 


7 If  I " (t)""]  p- -i  “ 


Proof;  We  rewrite  6v„  in  the  following  way: 


(A.  2) 


6Vn  = JJ  K(r,r,)  P,  (eos  jt)  d o 


where 


(A.  3) 


(r.r.)  = J k (r.rO  drj 


Thus  we  obtain: 


(A.  4) 


K (r,  r.)  = J ~ k (r,  rO  d ri  - k (r,  r) 
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where 
(A.  5) 
and 
(A.  6) 


k (r,r)  = 0 


By  combining  (A.  1)  through  (A.  6)  we  arrive  at  the  proposition. 


Proposition  2; 


Given: 


0 if  r s n 

K(r,r.)  = r { (n-1)  . if  r<r. 

L n+3  n-2 

irs/ri^“  1 

' * J + i Bn  (r,/r)  « if  r < Ts  , n = 2 


K (r,  rs)  can  be  e^qianded  into  the  following  series: 


(A.  7) 


where 


K (r,  r.)  = (2n+l)  H [^1  + + . . .] 


0 if  r a r. 

H = { 

'*  r.  - r otherwise 


Proof:  n / 2; 


(i)”^  1 . (««,  . lHl2p2  (f  )%  (a)%  . . . 


= 1 - ,.-2,  « (i)‘.  (l-J . ... 
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This  formula  satisfies  (A.  7)  for  n = 2. 


where 


implies 


M(M-H)  r IV-  Bi-3 
2M+1  L 2M-1 


IVt-3- 

2M+3 


where 


Proof:  We  have 


6^(N)  - J (to^-  t®)  y Pn  (t)  P„  (to)  d t 


From  Churchill  (1963,  p.  214),  we  obtain 


Thus  we  arrive  at: 


(to  + t)  (to)R<(t)  - Bii+i  (t)  (to)  d t 


6Ag(N)  = (N+1) 


[t  I\+i(to)B«(t)-toIVi(t)R<(to)]dt,  N = odd 


2 (N+1)  X 


[to  I\-.i  (to)  B-(t)  - 1 Pn*i  (t)I\  (to)l  d t , N = even 


In  the  last  derivation  we  have  used  the  fact  that  Bj  is  odd  if  N is  odd  and  even 
if  N is  even. 

We  use  the  following  substitution  (see  Churchill,  1963,  p.  206): 


(A.  11) 


t (t)  = 


1 

2N+1 


((N+1)  B(+i(t)  + N Bg-i(t)] 


Then  we  obtain  for  odd  N : 


(A.  12)  6Ag(N)  = 2 r tI\+i(to)PN-i(t)-I\+i(t)PN_i(to)ldt 

ZN+1  ^ 


PN4a  - Vti 


^ 2 (N+1)N  r PNj.Bi.-g. 

2N+1  L 2N-1  2N+3 


where  we  have  used  the  abbreviation: 


= Ptg  (to) 


(A.  12)  was  obtained  by  the  following  relations  from  Churchill  (1963,  p.  207): 


(A.  13) 
and 


(2n+l)  J P„  (t)  d t = P„i  (t)  - P„_i(t) 


Rk+i  (0)  = 0 

For  N even  we  obtain  in  the  same  way: 

^ ' 2N+3  L 2N+1  2N+5  J 
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IMittiiMMiagBiiBi 


implies 


where 


M(M+1) 

2M+1 


2M+3 


M 


N if  N is  odd 


N+1  if  N is  even 


and 


0 s t s to 


Proof:  From  the  derivation  of  Prop.  3,  we  obtain  (for  odd  N) 


(A.  14) 


where 


6V(N)  = J (to^-t^)  ^(t)  dt 


gN  = 2 ^1^  {I\^i(to)  R--i(t)-  R.«(t)R,-i(to)l 


Integrating  (A.  14)  by  parts,  we  obtain: 


6V(N)  = [(to®-  t®)  gNlo®  + 2 J t gN  (t)  d t 
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The  first  term  is  zero  because  gN  (0)  — 0 and  for  the  second  term  we  may 
use  the  mean  value  theorem  of  integral  calculus  (Bartle,  1964,  p.  303). 
The  result  is; 


(A.  15) 


where 


6V(N)  = 2 J t g«(t)  d t = g„(t) 


0^  ts  to 


In  the  same  way  the  proposition  is  proved  for  even  N. 


Corollary  4;  to  = 1 - 6v(N)  = 0 for  N ^ 4. 

Proof?  We  may  integrate  the  first  part  of  (A.  15)  one  more  time  by  parts. 
Then  we  obtain: 


where 


6 V(N)  = 2 to  Gn  (to)  - 2 J Gn  (t)  d t 


Gn  (t)  = J gN  (t)  d t 


From  the  expression  for  g^  in  Proposition  4 and  (A.  13)  it  follows  that; 


gN(l)  = Gn(1)  = J GN(t)  dt  = 0 forNi4 


which  implies 


6V  (N)  = 0 tor  to=  1 


A.  3 Estimates  of  the  Optimum  Degree  of  E^qpanslon  from  Surface  Mean  Anomalies 


r 


We  assume  that  we  have  a complete  coverage  of  mean  gravity 
[ anomalies  (Ag)  over  blocks  of  size  6°  x 6°  on  the  mean  earth  sphere.  As 

I we  have  a finite  set  of  mean  anomalies,  (bis  mean  gravity  field  can  always 

\ be  expanded  into  a finite  series  of  spherical  harmonics.  If  the  number  of 

[ coefficients  (k)  equals  the  number  of  mean  anomalies  (i)  there  is  a unique 

^ set  of  coefficients  An. , Bm  that  satisfies  all  equations: 


(A.  16) 


n=:0  BsrO 


(A„,  cosmX  + Bb,  sinmX)  Pb,  (sin<p)  d o 


where 


1 ~ 1,  2,  . . . , it 

N = degree  of  truncation 
Oi  = the  area  of  the  ith  block 


These  coefficients  can  therefore  be  determined  by  solving  a regular  matrix 
equation,  and  the  solution  will  give  the  exact  mean  value  over  each  block. 


Theoretically  we  are  not  limited  to  1=  k coefficients  in  order  to 
obtain  a soluticm  that  fits  all  Agi  . If  k > X (more  coefficients  than  observa- 
tions) the  solution  is  not  unique.  One  solution  is  thus  obtained  by  using 
condition  adjustment: 


YA  = 4^  - A = (YY^  )“^  Ag 


where 

A = vector  of  unknown  coefficients 

Y = vector  of  spherical  harmonics  (integrated  over  each  block) 

Ag  = vector  of  mean  anomalies 


Thus  there  are  different  sets  of  coefflcients  that  satisfy  the  observations  and 
there  is  no  evidence  that  one  set  is  superior  to  another. 
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In  a different  approach  we  eiqpand  Ag  Into  a complete  set  of  spherical 
harmonics: 


(A.  17) 


where 


^ = (So,  COS  m X + 


Ag  = ^ ^ (a„, 

n=0  BsO 


bga  sinmX)  Pb.  (sino) 


fa„,,  1 ..  cosmX.  _ 


If  this  series  is  truncated  it  will  not  be  complete  and  will  not  satisfy  the  observed 
mean  anomalies.  This  tendency  has  been  verified  by  Rapp  (1977,  p.  42).  Thus 
the  commonly  used  rule  of  the  optimum  degree  of  truncation: 


(A.  18) 


N - 180° 

^optimum  qo 


where  6°  is  the  block  size,  does  not  hold  for  mean  anomalies. 

We  may  then  ask  whether  (A.  18)  is  valid  for  estimating  point  anomalies 
(Ag).  We  assume  that  Ag  may  be  expanded  into  a convergent  series  at  each 
point  of  the  mean  earth  sphere: 


(A.  19) 


where 


Ag  = y Ag„ 


n = 0 


n 

Agn  = V (Cbi  cosmX  + Sn,  SinmX)  Pb,  (Sin  v>) 


1=0 


The  smoothed  field  (A.  17)  is  then  related  to  Ag  according  to: 

00 

(A. 20)  Ag  (X)  = JJ  B(x.y)  Ag(y)  d a,  = ^ AgB(x) 


B=0 


where  B (x*y)  is  the  integral  kernel  of  the  smoothing  operator  and  /S,  are  its 
eigen  values. 
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Assuming  that  each  block  has  the  same  size  as  a circular  disk  of 
spherical  radius: 


il>o  = (6  sin  6/n)^ 


we  obtain  from  Meissl  (1971): 


x«  y ^ cos  00 


B (x-rt  = { 2" 

0 otherwise 


1-COS0O 


(COS0o)-P 

2n+l 


Comparing  (A.  17),  (A.  19)  and  (A.  20)  we  obtain: 


(A. 21) 


{.3 = T.  {::} 


Thus  the  correct  spherical  harmonic  coefficients  of  ^ are  given  by  the 
coefficients  of  the  mean  anomalies  divided  by  • Theoretically  this 
estimate  of  A g is  improved  for  each  additional  coefficient  included  in  the 
series  (A.  19).  However,  the  estimates  of  the  coefficients  according  to 
(A.  21)  are  poor  for  higher  degrees,  because  )9„  approaches  zeit>.  Thus  we 
have  in  practice  to  limit  ourselves  to  a finite  number  of  coefficients. 

Let  us  now  assume  that  we  estimate  Ag  by  a truncated  form  of  (A.  20): 


fit 


The  error  of  this  estimate  is  then  given  by: 


” ^ m 

- Ag  = y Ag„  - 1)  - y Ag, 


As  the  spherical  harmonics  of  different  degrees  are  orthogonal  to  each 
other,  we  obtain  the  following  global  RMS  error; 


lliSg-  Agll  = [ya„^(l-  /3„)®  + f af  ]" 

»=0 

or 

00  N ^ 

{A.22)  iiAk-^ii=[y 

n=0  1^0 

where 

‘’'=  JJ 


We  define  Noptimum  as  the  value  of  N that  provides  the  minimum 
RMS  error  of  ^ [formula  (A.22)].  in  Figure  A.  1,  Noptimum  i®  given  as 
a function  of  the  block  size  6°  x 6°.  The  computations  are  based  on  the 
degree  variances  (a„^)  of  Tscheming  and  Rapp  (1974).  The  digram  clearly 
shows  that  the  optimum  degree  of  truncation  is  higher  than  is  indicated  by  i 

the  rule  (A.  18). 

Finally,  we  conclude  that  if  the  coefficients  of  (A.  16)  are  solved  from 
a system  of  linear  equations,  there  is  one  unique  solution  of  this  system  if 
the  number  of  equations  (£)  equals  the  number  of  coefficients  (k).  If  k > £ 
the  solution  is  not  unique.  We  may  not  conclude  that  one  set  of  coefficients 
is  generally  superior  to  another  set. 

If  the  spherical  harmonic  e;q)ansion  is  solved  by  integration,  the 
solution  will  converge  to  the  original  mean  anomalies  for  an  infinite  number 
of  coefficients.  Comparing  a truncated  form  of  this  series  with  point 
anomalies  we  arrive  at  the  conclusion  that  the  optimum  degpree  of  truncation 
is  higher  than  the  rule  of  thumb  (A.  18). 

i 


Figure  A.  1.  The  (^timum  degi^ee  of  truncation  of  a spherical  harmonic  expansion  of  the 
gravity  field  from  surface  mean  anomalies  of  block  size  0®  x 0®.  The 
computations  are  based  on  the  degree  variances  of  Tscheming  and  Rapp  (1974). 
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Table  A.  1.  Computed  error  coefficients  of  the  spherical  harmonic 
expansion  of  T at  sea  level.  A spherical  mean  earth 
is  anticipated.  The  coefficients  are  normalized: 

A = a.,  X r X 10® /GM , B = b„.  x r x 10®/GM 
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Table  A. 2.  Computed  error  coefficients  of  the  spherical  harmonic 
expansion  of  Ag  at  sea  level.  A spherical  mean  earth 
is  anticipated.  Unit;  mgal. 
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Table  A.  3.  Computed  error  coefficients  of  the  spherical  harmonic 
expansion  of  T for  r = 6371  km.  The  mean  earth  is 
anticipated  as  a level  ellipsoid  with  a = 6378. 140  km 
and  f = 298. 56.  The  coefficients  are  normalized; 

A = an,  X r X 10®/GM,  B = bn,  x r x 10®/GM. 
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